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Gravitational wave data analysis implications of TaylorEt inspiral approximants for 
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A new family of restricted post-Newtonian-accurate waveforms, termed TaylorEt approximants, 
was recently proposed for searching gravitational wave (GW) signals from inspiraling non-spinning 
compact binaries having arbitrary mass-ratios. One of the attractive features exhibited by these 
waveforms is that as their reactive post-Newtonian (PN) order is increased, their phase-evolution 
monotonically converges to that of waveforms produced by numerical relativity. The TaylorEt ap- 
proximant is different from the usual post-Newtonian ones, such as the TaylorTl, TaylorT4, and 
TaylorF2 approximants, in that it gives equal emphasis to both the conservative and the reactive 
parts of GW phase evolution. However, for the latter set of extensively employed PN-accurate inspi- 
ral templates the conservative phase evolution is somewhat dwarfed by the reactive part. We perform 
detailed fitting factor (FF) studies to probe if the TaylorEt (3.5PN) signals for non-spinning com- 
parable mass compact binaries can be effectually and faithfully searched with TaylorTl, TaylorT4, 
and TaylorF2 (3.5PN) templates in LIGO, Advanced LIGO, and Virgo interferometers. We observe 
that a good fraction of the templates, which by choice are from TaylorTl, TaylorT4, and TaylorF2 
(3.5PN) families, have FF ^ 0.97 and substantial biases for the estimated total-mass against the 
fiducial TaylorEt (3.5PN) signals for equal-mass systems. Both these observations can bear on the 
detectability of a signal. TaylorEt (3.5PN) signals with mass-ratios of a third or a quarter yield high 
FFs against those same template banks, but at the expense of inviting large systematic errors in 
the estimated values of their total mass and symmetric mass-ratio. In general, the aforementioned 
templates are found to be increasingly unfaithful with respect to a TaylorEt signal as one increases 
the total mass of the inspiraling system. We find that one way of improving the FF values is to 
allow the templates to have unphysical mass-ratios. However, this may result in a higher noise 
background and, therefore, reduce the detection confidence. We also observe that the amount of 
bias in the estimated mass varies with the (noise power spectral density of the) detector. This can 
be of some concern for multi-detector searches, which check for consistency in the estimated masses 
of concurrent triggers in their data. However, by modeling this variation it is possible to mitigate 
its effect on the detection efficiency. 

PACS numbers: 04.30.Db, 04.25.Nx 04.80.Nn, 97.60.Jd, 95.55Ym 



I. INTRODUCTION 



Stellar mass compact binaries, involving black holes 
and neutron stars, are the most promising sources of 
gravitational radiation for the operational and planned 
ground-based laser interferometric gravitational wave 
(GW) detectors. Gravitational wave signals from inspi- 
raling compact binaries are being searched in the de- 
tector data by using matched filtering [l| with several 
types of theoretically modeled inspiral templates [2, Q • 
A good resource for computing these templates in soft- 
ware, which is being actively used by the LIGO Scientific 
Collaboration (LSC) [J| and the Virgo Collaboration 5] 



'Electronic address: lsrikanta@wsu.edul 

t Currently at Max- Planck- Institut fiir Gravitationsphysik, Albert- 
Einstein- Institut, Callinstr. 38, 30167 Hannover, Germany 
-f Electronic address: a.gopakumar@uni-jena.de| 



§ Electronic address: |ni.tessmer@uni-jena.de| 



for analyzing detector data, is the LSC Algorithm Li- 
brary (LAL) 0. 

As such, the construction of these search templates re- 
quires two crucial inputs from the post-Newtonian (PN) 
approximation to general relativity, appropriate for de- 
scribing the dynamics of a compact binary during its in- 
spiral phase. These are the 3PN accurate dynamical (or- 
bital) energy £{x) and the 3.5PN accurate expression for 
the GW luminosity C{x) Q, both of which are usually 
expressed as PN series in the gauge invariant quantity 
X = {G moj / c^)'^/'^ , where m and w are the total mass 
and the orbital angular frequency of the binary, and the 
familiar symbols G and c denote the universal gravita- 
tional constant and the speed of light in vacuum, respec- 
tively. Recall that the 3PN accurate expression for £{x) 
provides corrections to the Newtonian orbital energy to 
the order of (w/c)®, where v is the orbital speed. Further, 
the currently employed search templates only require the 
Newtonian contributions to the amplitude of GW polar- 
izations, h^(t) and ^x (t). However, expressions for h-^-{t) 



and hx (t) that include the 3PN ampUtude corrections are 
available in Ref. [^ and are being used to develop ampli- 
tude corrected templates for GW inspiral searches. 

With the help of the two aforesaid PN-accurate inputs, 
one can construct two distinct classes of inspiral GW 
templates. The templates belonging to the first category 
require 
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and PN-accurate prescriptions for the reactive evolution 
of x{t). Such templates are usually referred to as adi- 
abatic inspiral templates, and all PN-accurate inspiral 
templates that LAL employs are of this type. In this pa- 
per, we consider from this class time-domain templates 
of the TaylorTl [^ and the Numerical Relativity (NR) 
inspired TaylorT4 |9|] families and frequency-domain tem- 
plates of the TaylorF2 family [l^. For all three families, 
we incorporate radiation reaction effects to the (relative) 
3.5PN order (see Eqs. (0), dJ), and ^ below). Due to 
their use of the x-based phase evolution expression in Eq. 
(HI), it may be argued that these templates model GWs 
from compact binaries inspiraling under PN-accurate ra- 
diation reaction along exact circular orbits [lOJ ] . 

A new class of inspiral approximants introduced in 
Rcf. 103, termed as TaylorEt, requires PN expansion 
for d(j}/dt in terms of the orbital binding Energy to de- 
rive the temporal GW phase evolution. This alternative 
phasing prescription models GWs from compact binaries 
inspiraling under PN-accurate reactive dynamics along 
PN-accurate circular orbits [lOj. In other words, the 
TaylorEt approximant explicitly incorporates the secular 
contributions to GW phase evolution appearing at the 
IPN, 2PN, and 3PN orders. Contrastingly, in the case 
of x-based adiabatic inspiral templates and due to the 
use of Eq. U), the above mentioned conservative (and 
secular) contributions to the GW phase evolution do not 
appear before the radiation reaction kicks in at the ab- 
solute 2.5PN order. It should be noted that the cost of 
computing TaylorEt templates is comparable to that of 
TaylorTl/T4 templates. 

In this paper, we study how effectively and faith- 
fully the TaylorTl, TaylorT4, and TaylorF2 inspiral tem- 
plates, at 3.5PN order, can capture a GW signal mod- 
eled using the TaylorEt approximant of the same order. 
The main motivation for using the latter as the fiducial 
signal originates from the observation that the TaylorEt 
approximant is an appropriate zero-eccentricity limit of 
GW phasing for compact binaries inspiraling along PN- 
accurate eccentric orbits [1^. We quantify our results 
by computing fitting factors (FF) following prescriptions 
detailed in Refs. [Il|, [I4I , and inherent systematic errors 
in the estimated value of m and the symmetric mass- 
ratio, rj, for the various search templates, relevant for 
the initial LIGO (heretofore referred to as "LIGO" ) , Ad- 
vanced LIGO (or "AdLIGO") and Virgo detectors. We 
conclude that it is desirable to incorporate the TaylorEt 
approximant at 3.5PN order into LAL to minimize pos- 



sible loss of inspiral events. Further, one might also view 
this work as an exercise in assessing the effects of us- 
ing inspiral templates from different representations on 
a GW signal's detectability and parameter estimation in 
earth-based detectors. Similar assessments of systematic 
errors on GW searches in LISA and Virgo were made by 
comparing inspiral templates of different PN orders from 
the same representation in Refs. 13] and [lj|, respec- 
tively. 

The plan of the paper is as follows. In the next sec- 
tion. Sec. ini we provide explicit PN-accurate equations 
required for constructing TaylorTl, Taylor T4, TaylorF2 
and TaylorEt templates having 3.5PN accurate reactive 
evolution. Section IIIII explains how we perform our FF 
computations and tabulates their values, along with the 
associated systematic errors in m and 77, for our different 
templates. We briefly discuss the implications of these 
results on the on-going searches in real interferometric 
data. We conclude in Sec. IIVI bv providing a brief sum- 
mary and future directions. 



II. PHASING FORMULAE FOR VARIOUS 
INSPIRAL TEMPLATES 



The PN approximation to general relativity is expected 
to describe accurately the adiabatic inspiral phase of a 
comparable mass compact binary Q . During this phase, 
the change in the orbital frequency over one orbit may 
be considered to be tiny compared to the mean orbital 
frequency itself. For compact binaries, having negligi- 
ble eccentricities, the adiabatic orbital phase evolution 
can be accurately described with the help of 3PN and 
3.5PN accurate expressions for the orbital energy and 
the GW luminosity, respectively, available in Refs. [3]. 
While employing a; as a PN expansion parameter, there 
exist several prescriptions to compute the adiabatic GW 
phase evolution. Each prescription, termed a PN approx- 
imant, provides a slightly different GW phase evolution 
and, correspondingly, a different inspiral template family. 
Following Ref. [g], we first list the equations describing 
the TaylorTl and the TaylorF2 approximants, which are 
regularly employed by various GW data analysis groups. 
The time-domain TaylorTl approximant is given by 
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where the proportionality constant in Eq. (j2ap may be 
set to unity for our analysis. To construct the TaylorTl 
3.5PN order adiabatic inspiral templates, one needs to 
use 3.5PN accurate C{x) and 3PN accurate S{x), respec- 
tively. The explicit expressions for these quantities, ex- 



tracted from Refs. [3], read 



C{x) 



''^'^\Ul 



5G 

44711 



1247 35 



336 12 



9271 



■V- 



65 



■77 



TT X 



9072 
583 
"24 
1712 

T05"^" ^ 7776 
775 o 1712 

■77^ 



■V 



X + Attx^/"^ 
8191 



504 ' 18 
5/2 . [6643739519 



69854400 



324 
214745 



/ 134543 41 

48 

In (4 V5) 



TT^ I 7; 



105 
193385 



672 

16^ 

3 

94403 
" 3024 
16285 

504 



E{x) 



1728 

„2 



■n 



,7/2 



rjmc 



24:^ 



x< 1- 



9 + r; 
35 



205 
"96" 



3024 
1 

12 
675 
"64" ~^ 5184 
34445 \ 



(3a) 



27 19 



155 
"96" 



576 



- . 



(3b) 



where 7 is the Euler constant and rj = ^/m, with /i being 
the reduced mass of the binary. 

The frequency-domain TaylorF2 approximant at 
3.5PN order, extracted from Ref. [l^, reads 
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where v — {Gtt m f/c^)^'^, and tc and (j)c are the fiducial 
time and phase of coalescence, respectively. The exphcit 
expressions for the PN coefficients ctfc are 
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where wiso is the speed at the last stable orbit, which we 
take to be at 6 Gm/c^. 

Recently, Ref. ^9] introduced another Taylor approxi- 
mant, termed TaylorT4. This approximant is obtained 
by Taylor expanding in x the right-hand side of Eq. (j2cp 
for dx/dt and truncating it at the appropriate reactive 
PN order. This approximant at 3.5PN order has an in- 
teresting (and accidental) property that was discovered 
due to the recent advances in Numerical Relativity (NR) 
involving coalescing binary black holes 16]. It was ob- 
served in Ref. 9] that the NR-based GW phase evolution 
for an equal-mass binary black hole agrees quite well with 
its counterpart in TaylorT4 approximant at 3.5PN order. 
Specifically, Ref. @] observed that the accumulated GW 
phase difference between TaylorT4 waveforms at 3.5PN 
order and NR waveforms agrees within 0.06 radians over 
30 wave cycles and matched at a; ^ 0.215. The time- 
domain TaylorT4 approximant at 3.5PN order is speci- 
fied by 
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It should be noted that the TaylorF2 waveform in Eqs. 
dU is the Fourier transform of h{t), given by Eqs. © 
above, computed with the help of the stationary phase 
approximation [l7|; we speculate that this is the reason 
that the TaylorT4 approximant is not directly employed 
in LAL (which already uses the TayIorF2 approximant). 
A close inspection of various time-domain adiabatic 
inspiral templates available in LAL reveals that they all 
invoke Eq. ^. These template families are different 
from one another only in the manner in which they in- 
corporate the reactive evolution of x{t). For example, 
PadeTl time-domain inspiral templates are constructed 
by invoking a specific Pade resummation for the right- 
hand side of Eq. (pc|) . Therefore, we may state that var- 
ious x-based inspiral templates provide slightly different 



GW phase evolution by perturbing a compact binary in 
an exact circular orbit, defined by Eq. ([1]), by different 
prescriptions for the reactive evolution of x{t). This is 
the main reason behind the observation that these tem- 
plates model GWs from compact binaries inspiraling un- 
der PN-accurate radiation reaction along exact circular 
orbits. 

Interestingly, it is possible to construct, in a gauge- 
invariant manner, inspiral GW search templates that do 
not require working in terms of the x variable. The 
TaylorEt approximant [lO| employs the orbital binding 
energy in lieu of that variable to describe PN-accurate 
adiabatic GW phase evolution. Hence, it requires an ap- 
propriate PN expansion for d(f>/dt in terms of the orbital 
binding energy. Accordingly, it can be argued that the 
TaylorEt approximant models GWs from compact bi- 
naries inspiraling under PN-accurate radiation reaction 
along PN-accurate circular orbits. 

The TaylorEt approximant at 3.5PN order is defined 
by 



h{t) (x£{t) cos2(/)(t) 

at Gm 

201 11 2 

■'^+128^ 



(7a) 



891 
128 



dm 

dt 



64 

205 
"64" 



2 [41445 /309715 
^ + mo. l^ 3072 



ttMt^ 



64 . r 



1215 2 
1024 ' 

.336 



1024 

45 

1024 



V' 



e 



(7b) 



|,)^ 



5 ^ ^„„„ 

/ 117857 12017 5 ^ 
^ \ 18144 2016 ^^2^ / 
,,/o [37999588601 

C' + 



inC 



:3/2 



e 



4913 
672 



177 



-V 



TTl 



^(4 7?) 



279417600 



1712 16 TT^ 



-7 



105 \ ^^J 105 

^369 2 24861497\ 488849 2 



^V32 

85 V 
-64^^ 



TT^ 



613373 2 
12096 ' 



72576 J ' 16128 
,„ [129817 3207739 
^ [ 2304 48384 ' 



(7c) 



where ^ = — 2£//ic^. A close inspection of Eqs. ([7]) re- 
veals that the above inspiral h{t) is obtained by perturb- 
ing a compact binary in a 3PN accurate circular orbit, 
defined by Eq. (|7bp. by radiation reaction effects at 3.5PN 
order, given by Eq. ((7c)l . The explicit use of PN-accurate 
expression for d(p/dt allows us to state that the TaylorEt 
approximants model GWs from compact binaries inspi- 
ralling under PN accurate reactive dynamics along PN 
accurate circular orbits. 

Importantly, a recent study of the accumulated phase 
difference between NR waveforms on the one hand and 
TaylorEt, TaylorTl, and TaylorT4 waveforms on the 



other hand reveals the following characteristics [l8|. In 
the interval x ~ 0.127 to a; ^ 0.215 this difference for the 
TaylorEt approximant at 3.5PN order is Scj) ^ —1.18 ra- 
dians, which is more than what is found for the TaylorTl 
and TaylorT4 counterparts (with Scj) ^ 0.6 and 0.06 ra- 
dians, respectively). However, significantly, TaylorEt is 
the only approximant studied so far that exhibits mono- 
tonic phase convergence with the NR waveforms when 
its reactive PN order is increased [l^. Recall that so- 
phisticated Pade approximations are required to make 
x-based Taylor approximants converge monotonically to 
the h(t) obtained from numerical relativity in the 77 = 
case [121 . In the context of the present paper, the analysis 
detailed in Ref. [l^l also suggests that TaylorEt approxi- 
mant at 3.5PN order remains fairly accurate in describing 
the inspiral h{t) even near the last stable orbit. These 
properties make it worthwhile to study the data analysis 
implications of TaylorEt approximants. 

Another motivation for using the TaylorEt approxi- 
mant to model the expected inspiral GW signal is as 
follows: With the help of Refs. [io|)llS|j it can be argued 
that TaylorEt is an appropriate approximant resulting 
from the zero-eccentricity limit of GW phasing of com- 
pact binaries inspiraling along PN-accurate eccentric or- 
bits. By contrast, the construction of the usual adiabatic 
inspiral templates requires redefining the right-hand side 
of Eq. (|7bp to be c'^x^/'^ /Gm, which can not be extended 
to yield GWs from precessing and inspiraling eccentric bi- 
naries as obtained in Ref. [1^. Therefore, the TaylorEt 
approximant can be expected to closely model GW sig- 
nals from inspiraling compact binaries, which realistically 
will not move along exactly circular orbits. The above 
statements are based on Ref. J20:| that, while restricting 
radiation reaction to the dominant quadrupole contribu- 
tions, demonstrated the undesirable consequences of re- 
defining the right-hand side of Eq. (|7b|) at 2PN order to 
be c^x^/"^ /Gm. 

Currently, for the low-mass binary signal searches 
the LSC usually employs templates based on Taylor Tn 
(where n=l, 2, and 3) and TaylorF2 approximants f6l.l27J. 
Therefore, it is important to probe if some of these tem- 
plates can capture inspiral signals modeled on the Tay- 
lorEt approximant. This is what we pursue in the next 
section. 



III. FITTING FACTORS 

Inspiraling compact binaries are the most promising 
sources of GWs for LIGO/Virgo. Detailed source popu- 
lation synthesis studies suggest that achieving an appre- 
ciable event rate, of at least a few compact binary coa- 
lescences per year, is possible if one could hear sources 
in the far reaches of our local super-cluster and beyond 
|2ll [23 |. Such an endeavor necessitates the ability to 
detect signals with relatively low signal-to-noise ratios 
(SNRs), even with second generation detectors, such as 
AdLIGO. Let the GW strain from a non-spinning com- 



pact binary be denoted by h{t; A), where A represents the 
signal parameters, namely, m, rj, tc, and (f>c, or an alter- 
native set of transformed coordinates in that parameter 
space. If a detector's strain-data is denoted by s{t) and 
its noise power-spectral-density (PSD) by S'„(/), then 
the SNR when filtering the data with template h{t; A') is 



SNR = 



,\h{X')) 



v/(MA')IMA')> ' 



(8) 



where A' symbolizes the template parameters, which 
need not be the same as the parameters of a signal em- 
bedded in the data, and the inner product {a\b) is defined 
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Above, a(/) is the Fourier transform of a{t) and the as- 
terisk denotes complex conjugation. Using Eq. ([8]), it can 
be shown that the quantity 
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also known as the "match" , is useful in describing how 
well two normalized waveforms, not necessarily from the 
same template family, overlap [23|. 

For the problem of detecting a GW inspiral signal, 
the prevailing sentiment in the community is that it is 
not as essential to search with a template bank that is 
an exact representation of the signal, as it is to search 
with an approximate one that can filter the data in real- 
time, provided its expected maximal match with a signal 
from anywhere in the parameter space is above a desired 
threshold. In other words, it should be possible to ob- 
tain a sufficiently large 'match' with a family of templates 
having A' ^ A. It is often stressed that this faithlessness 
of a template in accessing the signal parameters does not 
concern the detection problem per se, but that it affects 
the parameter-estimation problem, which can be tack- 
led a posteriori, i.e., after the transient signal has been 
detected and localized in time. The effectiveness of a 
template family, say, h{\'), in detecting the targ et signal 
g{X) is quantified by the fitting factor (FF) |ll| 



FF(A) 



max M (A, A') 

A' 



(11) 



If a template bank provides near-unity FF values for a 
given signal, it is considered to be effectual in detecting 
it [13. 

Employing an approximate template bank results in a 
drop in event rate by (1 — FF ) for a homogeneous dis- 
tribution of sources. This is easily seen when one realizes 
that the FF is a measure of the fractional loss in SNR 
(which scales inversely with source distance) stemming 
from using such a bank. So, e.g., a FF of 90% results 
in a 27% loss in event rate in any given detector. The 
expected rate in LIGO or Enhanced LIGO, which is a 



proposed upgrade in sensitivity of LIGO by roughly a 
factor of two while makin g m inimal changes to the shape 
of the LIGO noise PSD [2J], is very low, i.e., realisti- 
cally, less than one event in a few years. Therefore, a FF 
of 90% can potentially subvert a detection in the era of 
first-generation detectors. This is why a FF > 97% is so 
desirable. 

The values of the fitting factors for a couple of tem- 
plate banks against the TaylorEt 3.5PN waveforms are 
given in Tables [TTl lllli and IIVI The FFs were computed 
using two separate codes. One of these employs LAL [Q], 
which is used by the LSC in its inspiral searches, and the 
other is a home-grown code that extensively uses routines 
from Numerieal Reeipes 25]. Both codes have the ability 
to compute the FF in Eq. pT|) as well as the more con- 
servative (or lower) minimax match, detailed in Ref . [12| . 
The latter is obtained by minimizing the FF of Eq. pTJ) 
with respect to the coalescence phase of the target wave- 
form. The numbers presented in the tables are FFs (and 
not minimax matches) and, therefore, are larger than val- 
ues that are realistically achievable with the above listed 
inspiral template banks, available in LAL, and the Tay- 
lorT4 template bank. 

Importantly, the first few detections will likely require 
validation from more than one detector, which implies 
that in addition to being effectual a template bank also 
must be faithful [l2] ■ The latter requirement means that 
the parameter values of the best matched template are 
allowed to differ from (a subset of) those of the signal 
only by acceptably small biases. This is because unless 
these systematics are modeled for, the same signal can 
be picked up by templates with parameter values dif- 
ferent enough in two (or more) detectors so as to fail a 
parameter- value coincidence test |23|. We infer that dif- 
ferences in the estimated masses, illustrated in Tables |T] 
- IIVI between different comparable class detectors, such 
as AdLIGO and Virgo, are due to their different noise 
PSDs. 

Based on the tables and figures presented here, a few 
observations are in order. First, a good fraction of equal- 
mass compact binary templates, which are chosen to be 
from TaylorTl, TaylorT4 (presented only in the figures), 
and TaylorF2 (presented only in the tables) 3.5PN fam- 
ilies, have FF ^ 0.97. They also show substantial biases 
for the estimated total-mass against TaylorEt (3.5PN) 
signals as long as the symmetric mass-ratio of templates 
is limited to 77' < 0.25, which is the upper-limit for physi- 
cal signals. Note how in the first row of plots in Fig. [5]the 
FF first decreases as m is increased before recovering to 
higher values eventually. This behavior can be explained 
by the fact that in any given signal band the TaylorEt 
approximant has a greater number of GW cycles than the 
x-based templates of the same mass system. This means 
that a:-based templates with m' < m and rj' > rj are 
more likely to provide a higher match than the one with 
m' = m and 77' = 77. However, since for the equal-mass 
signals in Table [TTl and Fig. [2]we restrict the templates to 
have rf < 0.25, the templates yielding the highest match 



saturate this bound and attain rj' — r] — 0.25. This wave- 
based argument alone also implies that the highest match 
will decrease with increasing m. This is because while 
decreasing m' increases the number of template cycles, 
which helps in improving the match, decreasing it too 
much adversely affects the match by lowering the tem- 
plate /iso and, thereby, decreasing the integration band. 
The reason why the FF values eventually regain high val- 
ues at large m is that there the number of wave cycles is 
small and it is easier to obtain larger fits on signals with 
a smaller number of time- frequency bins. 

Compact binaries with mass-ratios smaller than unity 
can yield high FFs, but at the expense of introducing high 
systematic errors in estimating the values of m and rj. 
The high FFs can be explained by the fact that unlike in 
the case of equal-mass signals, here rj' can exceed 77, sim- 
ply because 77 < 0.25. For illustration purposes, we con- 
sider signals with two different values of the mass-ratio, 
namely, q = 7711/7712 = 1/3 and 1/4. In general, our in- 
spiral templates are found to be fairly unfaithful with re- 
spect to the fiducial signal for these cases. Moreover, Ta- 
bles [TT] and Fig. [T]show that the TaylorTl template with 
the maximum match almost always has rj' = 0.25. This 
arises from restricting the template banks to have physi- 
cal values of the symmetric mass-ratio, namely, rj' < 0.25. 
As shown in Fig. [TJ this can be seen from the fact that 
the match in Eq. (jlOp . also known as the ambiguity func- 
tion in this context, has a sharp wedge that rises with 
increasing rj' and attains its maximum beyond rj' = 0.25. 

The above observation prompted us to compute FFs 
in Tables IVl with 77' > 0.25. We find that it requires rj' 
to be as high as 0.35 for the FF to attain values at least 
equal to 97% for the total-mass range considered here 
(i.e., m < 40Mq), albeit, at the cost of large biases in the 
estimated values of both m' (up to almost 20%) and 77' 
(up to about 40%). In our opinion, while this allowance 
increases the match, it does not necessarily translate into 
an increase in the detection confidence. This is because 
an expanded range in rj' has the potential to increase the 
false-alarm rate. To test what the effective gain is it is 
imperative to include templates with 77' > 0.25 in signal 
simulation studies involving real interferometric data. 

As presented in Tables UlllIIIl the systematic errors 
also show variation with respect to the shape of the de- 
tector noise power spectral density. This implies, e.g., 
that the estimated value of the total mass of a signal in 
LIGO and Virgo can disagree and, consequently, fail a 
sufficiently stringent mass-consistency check in a multi- 
detector search [23, [2^. To wit, in the search for inspiral 
signals in LIGO data from its third and fourth science 
runs by the LSC, the estimated chirp-mass, Aic, h^^ the 
three LIGO detectors was allowed to differ by O.O2OAf0. 
A comparison of the estimated total-mass values in Ta- 
bles In] shows that this window needs to be relaxed if the 
search involved both LIGO and Virgo detectors and if 
TaylorEt was indeed a more appropriate representation 
of a GW signal. For instance. Table |TT] shows that for 
771 = 5.0 Mq the measured m' in LIGO and Virgo are ex- 



pected to differ by as much as 0.062Mq, which amounts 
to AAic — O.O26Af0 (where we assumed, conservatively, 
that all the error in Ale arose from the error inm'). This 
is larger than the allowed window and can, therefore, fail 
a multi-detector mass-consistency test. The biases shown 
in Tables HIl onlv get worse as the total mass of the signal 
is increased. Note that this exercise is meant to serve 
as a guide for sources of systematic effects and how to 
deal with them; it is not clear if a concurrent search with 
LIGO and Virgo detectors with LIGO and Virgo design 
sensitivity curves, respectively, (as used for Table ITT)) is 
likely. (Nevertheless, it is more likely that the shapes of 
the respective curves are maintained in the next set of 
science runs, such as with the planned Enhanced-LIGO 
design. This is because the FFs depend on the shape 
of these curves and not the overall scale). It is possible, 
however, to use our studies to model the variation of the 
estimated parameter bias in real detector data so that 
the windows can be scaled and shifted appropriately to 
mitigate the effect on detection efficiency. 



IV. CONCLUSIONS 

In this paper, we investigated the GW data analysis 
implications of the TaylorEt approximant at the 3.5PN 
order. We limited our attention to the case of GW sig- 
nals from non-spinning, comparable mass, compact bina- 
ries in LIGO, AdLIGO and Virgo interferometers. With 
the help of detailed fitting factor computations, we com- 
pared the performance of three x-based inspiral tem- 
plates, namely, TaylorTl, TaylorT4, and TaylorF2 at 
3.5PN order, in detecting a fiducial TaylorEt signal of 
the same PN order. For the equal-mass binaries, we gen- 
erally obtain FF ^ 0.97 when restricting the above tem- 
plates to physically allowed mass-ratios. In the case of 
unequal mass binaries, it is possible to obtain high FFs 
with the LAL inspiral templates. However, the templates 
that provide those high FFs have substantially different 
values of m and rj compared to those of the fiducial Tay- 
lorEt signals. In all cases, templates giving high FFs 
have lower values of the total-mass parameter compared 
to their associated TaylorEt signals. This is due to the 
fact that in a given GW frequency band the TaylorEt 
approximant always provides more accumulated GW cy- 
cles than the x-based templates. Further, the systematic 
errors in m and 77 parameters of TaylorF2 templates are 
substantially higher than the statistical errors in those 
parameters reported in Ref. T5] . These observations lead 
us to believe that the unfaithful nature of the x-based in- 
spiral templates vis a vis the TaylorEt approximant may 
adversely affect the chances of detecting GW inspiral sig- 
nals assuming that the latter waveforms more accurately 
represent such a signal. 

To summarize, the present study shows that it should 
be worthwhile to include the theoretically motivated Tay- 
lorEt templates, which has a number of attractive fea- 
tures as detailed in Ref. [10,] , in the search for inspiral GW 



signals from non-spinning compact binaries in the data 
of ground-based broadband detectors. Further, this work 
should also be useful in assessing the effects of systematic 
errors arising from employing inspiral templates from dif- 
ferent representations on a GW signal's detectability and 
parameter estimation with earth-based detectors. 

In the literature, there exist gauge-dependent prescrip- 
tions for constructing inspiral h(t) that give equal em- 
phasis to both the conservative and the reactive orbital 
phase evolution, such as the Effective One-Body (EOB) 
approach J23 and the Semi- Analytic Puncture Evolution 
(SAPE) [30|. Recall that the conservative Hamiltonian 
relevant for the EOB scheme is in Schwarzschild-type co- 
ordinates, while for SAPE it is in the Arnowitt-Deser- 
Misner gauge. By contrast, PN-accurate TaylorEt based 
h(t) is fully gauge- invariant. Therefore, we are pursuing 
a study, similar to the one presented here, of comparing 
the effectualness and faithfulness of EOB- and SAPE- 
based inspiral waveforms vis a vis TaylorEt waveforms. 
We are also extending the present analysis by including 
spin effects with the help of a generalized version of fidu- 
cial signals and templates detailed in Ref . [31| . 
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TABLE L The analytic fits to the one-sided noise power-spectral densities, S„(f), of LIGO, Virgo and AdLIGO employed in 
this paper. The expressions for Sn{f) are expressed in terms oi y = f/fo, where the "knee frequency" /o takes values of 150Hz, 
500Hz and 215Hz for the LIGO, Virgo and AdLIGO, respectively. Notice that Sn{f) rises sharply above the seismic cut-off fs- 
We drop an overall scale factor from each of the expressions of Sn{f) below as it does not affect our fitting factor studies. 



Detector 


/. 


S„{f) (up to an overall scale) 


LIGO 


40 Hz 


1 (4.49 j/)-^'^ + 0.16 y-*-52 + 0.52 + 0.32 y^ I 


Virgo 


20 Hz 


|(6.23t/)-^ + 2j/-i + l + y4 


AdLIGO 


20 Hz 


■^7,-*-" 59.-2 1 lll(l-H2+y''/2) 1 

\y ^y + (i-^yV2) J 



TABLE II: Values of the fitting factors and the associated systematic errors in the total mass parameter relevant for the LIGO, 
Virgo and AdLIGO detectors while employing the TaylorEt approximant at 3.5PN order to model the fiducial inspiral signal. 
The search templates, extracted from LAL, belong to the TaylorTl approximant at 3.5PN order in all cases except when the 
"Detector" is denoted as "AdLIGO-F2"; in the latter case the templates are given by the TaylorF2 approximant at 3.5PN 
order. The total mass and symmetric mass-ratio parameters of the signal and templates are denoted by {m,rf) and {m',r]'), 
respectively. We list m' of the search template that gives the largest match for a given signal. In all the systems below, 77' 
was found to be 0.250 and is, therefore, not listed. The fiducial GW signal and search templates are terminated when the 
instantaneous GW frequency reaches the value corresponding to the last stable orbit. The low FFs reported in certain cases 
may be attributed to the fact that we have not allowed -q' to go beyond its realistic bound of 0.25. The maximum overlap is 
always for values of m' < m. The numbers below are for the equal-mass case: q = mi/m2 = 1, i.e., r] = 0.25. It is interesting 
to note that for small m, Virgo has the worst FF. This owes its cause to a flatter noise PSD at low frequencies which weights 
the phase difference between the templates and the signal the most there. Also, the FF values of TaylorF2 templates are 
larger than those of TaylorTl (computed for comparison only for AdLIGO). This is consistent with the observation in Fig. [2] 
where the Taylor T4 FF values are also higher than those of TaylorTl. (Recall our speculation that TaylorF2 approximants are 
essentially stationary phase approximations of TaylorT4 approximants.) 



mi (Mq) = m2 (Mq) 


Detector 


FF 


m' (A/q) 


1.4 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.96 
0.98 
0.95 
0.93 


2.7995 
2.8000 
2.7985 
2.7999 


3.0 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.90 
0.94 
0.94 
0.87 


5.9924 
5.9952 
5.9841 
5.9950 


5.0 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.87 
0.92 
0.94 
0.84 


9.9365 
9.9404 
9.9117 
9.9740 


8.0 


AdLIGO-Tl 

AdLIGO-F2 

LIGO 

Virgo 


0.87 
0.92 
0.95 
0.86 


15.7649 
15.7624 
15.6663 
15.8770 


10.0 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.89 
0.90 
0.96 
0.89 


19.5718 
19.6786 
19.3438 
19.7728 


15.0 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.90 
0.93 
0.98 
0.92 


28.8432 
29.4024 
28.3201 
29.3476 


20.0 


AdLIGO-Tl 

AdLIGO-F2 

LIGO 

Virgo 


0.92 
0.89 
0.98 
0.93 


37.9325 
37.5663 
36.8851 
38.6133 
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TABLE III: Fitting factor values and inherent parameter biases relevant for the LIGO, Virgo and AdLIGO interferometers, 
for several compact binaries with mass-ratio q — 1/3 (i.e., 77 = 0.1875)). The other details are as in Table HIl We observe that 
high FFs are always for templates characterized by m' < m and 77' ~ 0.25. It is clear that the best matched templates are 
highly biased with respect to their fiducial signals. 



mi {Mq) - m2 (Mq) 


m (Mq) 


Detector 


FF 


m' (Mq) 


v' 


3-9 


12 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.98 
0.97 
0.98 
0.97 


10.1630 
10.2269 
10.1588 
10.1596 


0.250 
0.248 
0.250 
0.250 


4-12 


16 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.97 
0.98 
0.98 
0.97 


13.5519 
13.5704 
13.5078 
13.5522 


0.250 
0.250 
0.250 
0.250 


5-15 


20 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.96 
0.97 
0.98 
0.97 


16.9127 
16.9609 
16.8320 
16.9299 


0.250 
0.250 
0.250 
0.250 


7-21 


28 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.96 
0.96 
0.99 
0.97 


23.6070 
23.7107 
23.4532 
23.6543 


0.250 
0.250 
0.247 
0.250 


10-30 


40 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.98 
0.95 
0.98 
0.98 


33.5956 
33.8244 
35.4216 
33.6106 


0.250 
0.249 
0.208 
0.250 



TABLE IV: Fitting factor values, relevant for the LIGO, Virgo, and AdLIGO, for several compact binaries having mass ratio 
q = 1/4 or 77 = 0.16. The other details are as in Table Ull and our conclusions are also similar to those reported in Table HITl 
Additionally, it is interesting to observe that across both Table HITl and the current table, as the total mass of the system is 
increased in LIGO, the bias in 77' there appears to decrease noticably. This is because LIGO has a higher fs than the other 
detectors considered here (see Table [I} . Correspondingly, the number of in-band cycles is considerably reduced in it for higher 
mass, to the extent that tuning only one parameter, m' , in LIGO suffices to attain higher FFs. 



mi (Mq) - 7772 (Mq) 


m (Mq) 


Detector 


FF 


m' (Mq) 


V 


3-12 


15 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.97 
0.97 
0.99 
0.95 


11.7544 
12.2512 
11.9641 
11.7448 


0.246 
0.228 
0.238 
0.246 


4-16 


20 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.98 
0.97 
0.99 
0.98 


15.5198 
16.0009 
15.5188 
15.5213 


0.250 
0.237 
0.250 
0.250 


5-20 


25 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.98 
0.96 
0.99 
0.99 


19.4193 
20.2636 
20.3615 
19.4180 


0.250 
0.231 
0.224 
0.250 


28-7 


35 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.96 
0.96 
0.98 
0.98 


27.2169 
27.6552 
28.7005 
27.2233 


0.250 
0.243 
0.217 
0.250 


8-32 


40 


AdLIGO-Tl 
AdLIGO-F2 

LIGO 

Virgo 


0.98 
0.95 
0.98 
0.98 


32.0353 
32.0962 
35.9969 
31.9093 


0.234 
0.235 
0.167 
0.236 
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TABLE V: The recomputed FF values for AdLIGO for compact binaries studied in Table [IT] Unlike in that table, here we 
allowed the symmetric mass-ratio parameter of the templates to vary beyond the physically allowed range, and up to 77' = 0.35. 
As shown below, doing so improves the FFs to be higher than or equal to the desirable lower-limit of 97%. However, this gain 
comes at the cost of high biases in estimated masses and, possibly, an increased noise contribution to the detection statistic. 



mi (Mq) 


FF 


m' (Mq) 


v' 


1.4 


0.99 


2.7350 


0.260 


3 


0.99 


5.5939 


0.288 


5 


0.99 


8.8575 


0.308 


8 


0.99 


13.4103 


0.341 


10 


0.98 


16.5826 


0.348 


15 


0.97 


24.7678 


0.350 


20 


0.97 


32.6322 


0.350 



12 




il'?il ^ " m'/m 



q=0.25 




FIG. 1: The plots display the match M (see Eq. llOp [12] as a function of the parameters (m', t)') of the search templates. The 
plot axes are r)' /r) and m'/m, where {m,,r]) are the values of the TaylorEt signal parameters. For the equal-mass case {q — 1), 
the match keeps rising as rj' is increased and attains its global maximum in the unphysical region rj' > 0.25. In the q = 1/4 
plot, the maximum value of M is around rj' ~ 0.25, which is much larger than the signal parameter value of r; = 0.16. The 
maximum match remains fairly high, at more than 0.95, on certain crests of the ridges. 
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FIG. 2: A collection of plots that summarizes a set of results from our fitting-factor studies involving the TaylorEt, Tay- 
lorTl/TaylorT4 approximants at 3.5PN order for equal-mass binaries. The first row provides plots of FF against the total 
mass of the inspiral signals and the second row shows fractional error (in percentage) in the estimated total-mass values as a 
function of the signal's total-mass. The thick and dashed lines denote results for the TaylorTl and TaylorT4 approximants, 
respectively. For high mass binaries, we observe high biases in the m' values. Above, we do not provide 77' versus m plots as 
the best matched filters always have rj' ~ 0.25. 
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FIG. 3: A set of plots summarizing our fitting-factor and parameter estimation results for the TaylorEt, TaylorTl approximants 
at 3.5PN order for compact binaries having q — 1/4. The first two rows show plots analogous to those presented in Fig. [2] 
The third row shows plots of r/' versus m. (Note that r; = 0.16 in all these plots.) Here too the thick and dashed lines are for 
the TaylorTl and TaylorT4 approximants, respectively. We clearly observe substantially higher biases for 77' that for m' values 
and in most cases, the best matched templates always have 77' ~ 0.25. 



